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Abstract 

The gauge invariant generation of an effective gluon mass proceeds through the weh-known 
Schwinger mechanism, whose key dynamical ingredient is the nonperturbative formation of lon- 
gitudinally coupled massless bound-state excitations. These excitations introduce poles in the 
vertices of the theory, in such a way as to maintain the Slavnov- Taylor identities intact in the pres- 
ence of massive gluon propagators. In the present work we first focus on the modifications induced 
to the nonperturbative three-gluon vertex by the inclusion of massless two-gluon bound-states into 
the kernels appearing in its skeleton-expansion. Certain general relations between the basic build- 
ing blocks of these bound-states and the gluon mass are then obtained from the Slavnov- Taylor 
identities and the Schwinger-Dyson equation governing the gluon propagator. The homogeneous 
Bethe-Salpeter equation determining the wave-function of the aforementioned bound state is then 
derived, under certain simplifying assumptions. It is then shown, through a detailed analytical and 
numerical study, that this equation admits non-trivial solutions, indicating that the QCD dynam- 
ics support indeed the formation of such massless bound states. These solutions are subsequently 
used, in conjunction with the aforementioned relations, to determine the momentum-dependence 
of the dynamical gluon mass. Finally, further possibilities and open questions are briefly discussed. 

PACS numbers: 12.38.Lg, 12.38.Aw, 12.38.Gc 
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I. INTRODUCTION 



The numerous large-volume lattice simulations carried out in recent years have firmly 
established that, in the Landau gauge, the gluon propagator and theghost dressing function 
of pure Yang-Mills theories are infrared finite, both in SU{2) [JQ and in SU{3) UQ. 
Perhaps the most physical way of explaining the observed finiteness of these quantities is 



the generation of a non-perturbative, momentum-dependent gluon mass [10|-[l5|, which acts 
as a natural infrared cutoff. In this picture the fundamental Lagrangian of the Yang-Mills 
theory (or that of QCD) remains unaltered, and the generation of the gluon mass takes 



16 



22| without violating 



place dynamically, through the well-known Schwinger mechanism 
any of the underlying symmetries (for further studies and alternative approaches, see, e.g., 
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H29|). 



The way how the Schwinger mechanism generates a mass for the gauge boson (gluon) 
can be seen most directly at the level of its inverse propagator, A~-'^(g^) = ^^[1 + in(g^)], 
where n(g) is the dimensionless vacuum polarization. According to Schwinger's funda- 
mental observation, if n(g^) develops a pole at zero momentum transfer (g^ = 0), then 
the vector meson acquires a mass, even if the gauge symmetry forbids a mass term at the 
level of the fundamental Lagrangian. Indeed, if n(g^) = m? /q^, then (in Euclidean space) 
A^^(g^) = + m? , and so the vector meson becomes massive, A~^(0) = m^, even though 



it is massless in the absence of interactions ((? = 0, 11 = 0) [18|, Il9l |. 

The key assumption when invoking the Schwinger mechanism in Yang-Mills theories, such 
as QCD, is that the required poles may be produced due to purely dynamical reasons; specif- 
ically, one assumes that, for sufficiently strong binding, the mass of the appropriate bound 
state may be reduced to zero 
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g|-|22l]. In addition to triggering the Schwinger mechanism, 
these massless composite excitations are crucial for preserving gauge invariance. Specifically, 
the presence of massless poles in the off-shell interaction vertices guarantees that the Ward 
identities (Wis) and Slavnov Taylor identities (STIs) of the theory maintain exactly the same 
form before and after mass generation (i.e. when the the massless propagators appearing in 



them are replaced by massive ones) 
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22| . Thus, these excitations act like dynam- 



ical Nambu-Goldstone scalars, displaying, in fact, all their typical characteristics, such as 
masslessness, compositeness, and longitudinal coupling; note, however, that they differ from 
Nambu-Goldstone bosons as far as their origin is concerned, since they are not associated 
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with the spontaneous breaking of any global symmetry lOj. Finally, every such Goldstone- 
like scalar, "absorbed" by a gluon in order to acquire a mass, is expected to act ually cancel 



out of the S"- matrix against other massless poles or due to current conservation |18H22|. 

The main purpose of the present article is to scrutinize the central assumption of the 
dynamical scenario outlined above, namely the possibility of actual formation of such mass- 
less excitations. The question we want to address is whether the non-perturbative Yang 
Mills dynamics are indeed compatible with the generation of such a special bound-state. 
In particular, as has already been explained in previous works, the entire mechanism of 
gluon mass generation hinges on the appearance of massless poles inside the nonperturba- 
tive three-gluon vertex, which enters in the Schwinger Dyson equation (SDE) governing the 
gluon propagator. These poles correspond to the propagator of the scalar massless excita- 
tion, and interact with a pair of gluons through a very characteristic proper vertex, which, of 
course, must be non vanishing, or else the entire construction collapses. The way to establish 
the existence of this latter vertex is through the study of the homogeneous Bethe-Salpeter 
equation (BSE) that it satisfies, and look for non-trivial solutions, subject to the numerous 
stringent constraints imposed by gauge invariance. 

This particular methodology has been adopted in various early contributions on this 
subject; however, only asymptotic solutions to the corresponding equations have been con- 
sidered. The detailed numerical study presented here demonstrates that, under certain 
simplifying assumptions for the structure of its kernel, the aforementioned integral equation 
has indeed non-trivial solutions, valid for the entire range of physical momenta. This result, 
although approximate and not fully conclusive, furnishes additional support in favor of the 
concrete mass generation mechanism described earlier. 

The article is organized as follows. In Section |TT] we set up the general theoretical frame- 
work related to the gauge-invariant generation of a gluon mass; in particular, we outline how 
the vertices of the theory must be modified, through the inclusion of longitudinally coupled 
massless poles, in order to maintain the Wis and STIs of the theory intact. In Section IIIII 
we take a detailed look into the structure of the non-perturbative vertex that contains the 
required massless poles, and study its main dynamical building blocks, and in particular 
the transition amplitude between a gluon and a massless excitation and the proper vertex 
function (bound-state wave function), controlling the interaction of the massless excitation 
with two gluons. In addition, we derive an exact relation between these two quantities and 
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the first derivative of tlie (momentum-dependent) gluon mass. Tlien, we derive a simple 
formula that, at zero momentum transfer, relates the aforementioned transition amplitude 
to the gluon mass. In the next two sections we turn to the central question of this work, 
namely the dynamical realization of the massless excitation within the Yang-Mills theory. 
Specifically, in Section |IV] we derive the BSE that the proper vertex function satisfies, and 
implement a number of simplifying assumptions. Then, in Section |V] we demonstrate through 
a detailed numerical study that the resulting homogeneous integral equation admits indeed 
non-trivial solutions, thus corroborating the existence of the required bound-state excita- 
tions. In Section I VI I we demonstrate with a specific example the general mechanism that 
leads to the decoupling of all massless poles from the physical (on-shell) amplitude. Finally, 
in Section IVIII we discuss our results and present our conclusions. 

II. GENERAL CONSIDERATIONS 

In this section, after establishing the necessary notation, we briefly review why the dy- 
namical generation of a mass is inextricably connected to the presence of a special vertex, 
which exactly compensates for the appearance of massive instead of massless propagators in 
the corresponding Wis and STIs. 

The full gluon propagator A^^(g) = 5'^^A^,^(g) in the Landau gauge is defined as 

A^,{q) = -tP^,{q)A{q') , (2.1) 

where 

P,u{q) = 9,u - , (2.2) 

is the usual transverse projector, and the scalar cofactor A(g^) is related to the (all-order) 
gluon self-energy n^,y(g) = P^u{q)Il{q'^) through 

A-\q')=q' + tn{q'). (2.3) 

One may define the dimensionless vacuum polarization n(g^) by setting n(g^) = g^n(g^) 
so that (12. 3p becomes 

A-\q') = q'[l + tn{q')] . (2.4) 

As explained in the Introduction, if 11 (g^) develops at zero momentum transfer a pole with 
positive residue m^, then A^^(O) = m^, and the gluon is endowed with an effective mass. 
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FIG. 1: The "one-loop dressed" gluon contribution to the PT-BFM gluon self-energy. White 
(black) circles denote fully dressed propagators (vertices); a gray circle attached to the external 
legs indicates that they are background gluons. Within the PT-BFM framework these two diagrams 
constitute a transverse subset of the full gluon SDE. 



Alternatively, one may define the gluon dressing function J{q^) as 



A-i(g^) = gV(g^) 



(2.5) 



In the presence of a dynamically generated mass, the natural form of A ^{q^) is given by 
(Euclidean space) 

^-\q') = q^J{q')+m\q'), (2.6) 

where the first term corresponds to the "kinetic term", or "wave function" contribution, 
whereas the second is the (positive-definite) momentum-dependent mass. If one insist on 
maintaining the form of (12. Sp by explicitly factoring out a g^, then 

m^(g^) 



A-^(g^ 



W + 



(2.7) 



and the presence of the pole, with residue given by m^(0), becomes manifest. 

Of course, in order to obtain the full dynamics, such as, for example, the momentum- 
dependence of the dynamical mass, one must turn eventually to the SDE that governs the 
corresponding gauge-boson self-energy (see Fig. [1]). In what follows we will work within the 
specific framework provided by the synthesis of the pinch technique (PT) 10|, l30N34| with 



the background field method (BFM) 35||. One of the main advantages of the "PT-BFM" 

f-energy n^j^(g), namely 



formalism is that the crucial transversality property of the gluon se^ 
q^Ii^jji^q) = , is maintained at the level of the truncated SDEs 12 



36| 



The Schwinger mechanism is integrated into the SDE of the gluon propagator through 
the form of the three-gluon vertex. In particular, as has been emphasized in some of the 
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literature cited above (e.g., 15j), a crucial condition for the realization of the gluon mass 
generation scenario is the existence of a special vertex, to be denoted by Va^i,{q,r,p) which 
must be completely longitudinally coupled, i.e. must satisfy 



P"'"(g)P'^''^(r)P-'-(p)KM.(g,r,p) = 



(2.8) 



We will refer to this special vertex as the "pole vertex" or simply "the vertex V" . 

The role of the vertex Va^u{(l-,^-,p) is indispensable for maintaining gauge invariance, 
given that the massless poles that it must contain in order to trigger the Schwinger mecha- 
nism, act, at the same time, as composite, longitudinally coupled Nambu-Goldstone bosons. 
Specifically, in order to preserve the gauge-invariance of the theory in the presence of masses, 
the vertex V^^i,(g, r, p) must be added to the conventional (fully-dressed) three-gluon vertex 
rQ,^j,(g, r, p), giving rise to the new full vertex, r^^j^(g, r, p), defined as 



(2.9) 



Gauge-invariance remains intact because F' satisfies the same STIs as F before, but now 
replacing the gluon propagators appearing on their rhs by a massive ones; schematically, 
A^^ — !■ A~^, where the former denotes the propagator given in (12.51) . while the latter that 
of (12. 6p . In particular, in the PT-BFM framework that we employ, the vertex F connects a 
background gluon (B) with two quantum gluons (Q), and is often referred to as the "BQQ" 
vertex. This vertex satisfies a (ghost-free) WI when contracted with the momentum of 
the background gluon, whereas it satisfies a STI when contracted with the momentum or 
Pu of the quantum gluons. In particular, 

g"r„^,(g,r,p) = p^J{p^)P^,{p)-r^J{r^)P,,{r), 

r^r,^,(g,r,p) = F{r') [q^J{q^)P^{q)H^,{q,r,p) - p^J{p')Pi:{p)H^^{p,r,q) 
p''T^,,{q,r-p) = F{p') [rV(r2)p;(r)^.4r,p, g) - J(g2)p^.(^)^^^(^^p^ ^ (2.10) 

where P(g^) is the "ghost dressing function" , defined as P(g^) = q^D{q^), Hy^ is the standard 
gluon-ghost kernel, and H is the same as H but with the external quantum gluon replaced 
by a background gluon. Similarly, J is the dressing function of the self-energy connecting a 



background with a quantum gluon; J is related to J{q^) through the identity |37|, [38 1 



J(g^) = [1 + G(g^)] J(g^) . 



(2.11) 



The function G{q'^) is the scalar co-factor of the gf^j, component of the special two-point 
function A^^{q), defined as 

Ku{q) = -tg^CA [ A;{k)D{q - k)H,^{-q,q - k,k) 
Jk 



Note finally that, in t 



he 



(all-order) relation 39l-l42 1 



= 9,MQ') + ^-^Liq')- (2.12) 
;andau gauge, G^q"^) and L^q"^) are linked to -F(g^) by the exact 



F-\q') = 1 + G{q') + L{q') , (2.13) 

to be employed in Subsection D. 

Returning to the nonperturbative vertex V, gauge invariance requires that it must satisfy 
the WI and STI of (12.101) . with the replacement k'^J{k) — )■ —m'^{k), e.g., 

g"V;^,(g,r,p) = m\r^)P^,ir) ~ m\p'')P^,{p) ; (2.14) 

exactly analogous expressions will hold for the STIs satisfied when contracting with the 
momenta r or p. Indeed, under this assumption, the full vertex F' will satisfy the same 
WI and STIs as the vertex F before the introduction of any masses, but now with the 
replacement q^J{q^) — ?■ q^J{q^) + ra^{q^). Specifically, combining the first relation in (I2.10p 
with f l2J4D . one obtains for the WI of T\ 

g"r;^,(g,r,p) = [r(g, r,p) + l^(g, r,p)]^^, 

= [pV(p2) - m\p')]P^,{p) - [rV(r2) - m\r^)]P^,{r) 

= A~\p')P,,ip) - A-^\r^)P,,{r) , (2.15) 

which is indeed the first identity in Eq. ( ]2.10p . with the aforementioned replacement 
A^^ —7- A^^ enforced. The remaining two STIs are realized in exactly the same fashion. 

It must be clear at this point that the longitudinal nature of Vap,v, combined with the 
WI and STIs that it must satisfy, lead inevitably to the appearance of a massless pole, as 
required by the Schwinger mechanism. For example, focusing only on the g-channel, the 
simplest toy Ansatz for the vertex is 

Vo.^^M.r.p) = '^[m\r^)P^,{r) - m\p')P^,{p)] , (2.16) 

which has a pole in q^ and satisfies (12.141) . Of course, poles associated to the other channels 
(r and p) will also appear, given that Vafj,u{q,r,p) must also satisfy the corresponding STIs 
with respect to and p" . 



III. THE POLE VERTEX: STRUCTURE AND PROPERTIES 



In this section we have a detailed look at the structure of the special vertex V. In 
particular, we identify the diagrammatic origin and field-theoretic nature of the various 
quantities contributing to it, and specify the way it enters into the SDE of the full vertex 
r', defined in Eq. (12. 9p . In addition, we will derive an exact relation between the most 
important component of this vertex and the derivative of the momentum-dependent gluon 
mass. 

A. General structure of the vertex V 

The main characteristic of the vertex V, which sharply differentiates it from ordinary 
vertex contributions, is that it contains massless poles, originating from the contributions of 
bound-state excitations. Specifically, all terms of the vertex V are proportional to l/g^, 

and products thereof. Such dynamically generated poles are to be clearly distinguished 
from poles related to ordinary massless propagators, associated with elementary fields in the 
original Lagrangian. 



a 




[ai] (0.2) (03) (04) (as) 



FIG. 2: The SDE for the BQQ vertex which connects a background gluon (B) with two quantum 
gluons (Q). 

In general, when setting up the usual SDE for any vertex (see, for example. Fig. [2]), a 
particular field (leg) is singled out, and is connected to the various multiparticle kernels 
through all elementary vertices of the theory involving this field (leg). The remaining legs 
enter into the various diagrams through the aforementioned multiparticle kernels (black 
circles in graphs 02-05 in Fig. [2]), or, in terms of the standard skeleton expansion, through 
fully-dressed vertices (instead of tree-level ones). For the case of the BQQ vertex that we 
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FIG. 3: (A) A diagram that does not belong to the standard kernel. (B) The gray kernel (regular 
part with respect to q, and the composite massless excitation in the g-channel. (C) The R part of 
the vertex. 

consider here [shown in Fig. |2] , it is convenient (but not obhgatory) to identify as the special 
leg the background gluon, carrying momentum q. Now, with the Schwinger mechanism 
turned off, the various multiparticle kernels appearing in the SDK for the BQQ vertex have 
a complicated skeleton expansion (not shown here), but their common characteristic is that 
they are one-particle-irreducible with respect to cuts in the direction of the momentum g; 
thus, a diagram such as the (a) of Fig. [3] is explicitly excluded from the (gray) four-gluon 
kernel, and the same is true for all other kernels. 

When the Schwinger mechanism is turned on, the structure of the kernels is modified by 
the presence of composite massless excitation, described by a propagator of the type i/q^, 
as shown in Fig. [3l The sum of such dynamical terms, coming from all multiparticle kernels, 
shown in Fig. H] constitutes a characteristic part of the vertex V, to be denoted by U in 
Eq. (13. 4p . namely the part that contains at least a massless propagator i/q^. The remaining 
parts of the vertex V , to be denoted by R in Eq. (13.51) . contain massless excitations in the 
other two channels, namely and l/p^ (but no l/q"^), and originate from graphs such 
as (C2) of Fig. [HI Indeed, note that the kernel (62) is composed by an infinite number of 
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diagrams, such as (ci), containing the full vertex JT'; these graphs, in turn, will furnish terms 
proportional to and [e.g., graph (02)]. 

In order to study further the structure and properties of the vertex V, let us first define 
the full vertex V^™"(g, r, p), given by 

yZ:{q,r,p) = gn^V^,M,r,p) , (3.1) 

with Va^u{(l,r,p) satisfying Eq. (12. 8p . Using a general Lorentz basis, we have the following 
expansion for Va^y{q,r,p) in terms of scalar form factors, 

Va^„.{q, r,p) = Viq^g^u + Vag^g^g,. + V^q^p^^Py + Viq^r^^qy + Vsg^r^p,, 

+ Var^^gc^y + Vjrarf^ry + Vgr^r^p^ + V^pyg^^, + Viop^p^^Pu ■ (3.2) 

According to the arguments presented above, Va^y{q,r,p) may be decomposed into 

K«/.i.(g,r,p) = Uaf,y{q,r,p) + Raf,u{q,r,p) , (3.3) 

with 

UaiMuiq, r, p) = qa (vig^„ + V2q^q„ + Vsp^p,, + V^r^q^ + V^r^p^ . (3.4) 

and 

Raf,u{q, r, p) = (vegau + Vrr^ru + y r„p^^ + r^r^ + Ygg^^^ + VioPaP,^ Pu ■ (3.5) 

All form-factors of U (namely V1-V5) must contain a pole 1/g^, while some of them may 
contain, in addition, 1/r^ and l/p"^ poles. On the other hand, none of the form-factors of R 
(namely Ve-Vio) contains 1/g^ poles, but only 1/r^ and l/p"^ poles. 

In what follows we will focus on Uanuiq, i^,p), which contains the explicit g-channel mass- 
less excitation, since this is the relevant channel in the SDE of the gluon propagator, where 
Va^i,{q,r,p) will be eventually inserted [graph (ai) in Fig. [T]. In fact, with the two internal 
gluons of diagram (oi) in the Landau gauge, we have that 

P'^'^(r)P^''^(p)K^,(g,r,p) = P^'^(r)P^'^(p)f/„^,(g, r,p) 

= P^'^{r)P'''%p)q4V,{q,r,p)g^, + V2{q,r,p)q,q,], (3.6) 

so that the only relevant form factors are Vi and V2. 
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FIG. 4: (A) The vertex U is composed of three main ingredients: the transition amplitude, 
la, which mixes the gluon with a massless excitation, the propagator of the massless excitation, 
and the (massless excitation)-(gluon)-(gluon) vertex. (B) The Feynman rules (with color factors 
included) for (i) the propagator of the massless excitation and (ii) the "proper vertex function", 
or, "bound-state wave function", B^i,. (C) The various -B{ } appearing in Eq. ()3.14p . 

At this point we can make the nonperturbative pole manifest, and cast Ua^iviq.r.p) in 
the form of Fig. HJ by setting 

Ua^^uiq,r,p) = Ia{q) (^-^^ B^^{q,r,p) , (3.7) 

where the nonperturbative quantity 

Bi„.iq, r, p) = Big^^ + B2qf,q^ + B^p^p^ + B^r^q^ + B^^r^^p^ , (3.8) 

is the effective vertex (or "proper vertex function" |l^) describing the interaction between 
the massless excitation and two gluons. In the standard language used in bound-state 
physics, B^^{q,r,p) represents the "bound-state wave function" (or "BS wave function") of 
the two-gluon bound-state shown in (63) of Fig. [31 as we will see in Section [TV] B^^ satisfies 
a (homogeneous) BSE. In addition, i/q"^ is the propagator of the scalar massless excitation. 
Finally, Ia{q) is the (nonperturbative) transition amplitude introduced in Fig. HI allowing 
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the mixing between a gluon and the massless excitation. Note that this latter function is 
universal, in the sense that it enters not only in the pole part V associated with the three- 
gluon vertex, but rather in all possible such pole parts associated with all other vertices, 
such as the four-gluon vertex, the gluon-ghost-ghost vertex, etc (see panel C in Fig. 
Evidently, by Lorentz invariance, 

Uq) = QJiq) , (3.9) 
and the scalar cofactor, to be referred to as the "transition function" , is simply given by 

/(«) = ^. (3.10) 

SO that 

V,{q,r,p) = I{q)(^-^^B,{q,r,p); j = l,...,5. (3.11) 

Note that, due to Bose symmetry (already at the level of V) with respect to the inter- 
change fi ■H- u and p ^ r, we must have 

Bi,2iq,r,p) = -Bi^2iq,P,r) , (3.12) 

which implies that 

5i,2(0,-p,p) = 0. (3.13) 

Finally, in principle, all other elementary vertices of the theory may also develop pole 
parts, which will play a role completely analogous to that of Vafj,u in maintaining the cor- 
responding STIs in the presence of a gluon mass. Specifically, in the absence of quarks, 
the remaining vertices are the gluon-ghost-ghost vertex, Fq, the four-gluon vertex JTa^^p, 
and the gluon-gluon-ghost-ghost vertex r^^, which is particular to the PT-BFM formula- 



tion. The parts of their pole vertices containing the 1/g , denoted by Ua, Ua^iup, and f/, 



respectively, will all assume the common form 



I 



f/.{...} = /«(^)5{...}, (3.14) 



where the various are shown in panel C of Fig. HI 



B. An exact relation 



The WI of Eq (12.141) furnishes an exact relation between the dynamical gluon mass, 
the transition amplitude at zero momentum transfer, and the form factor Bi. Specifically, 
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contracting both sides of the WI with two transverse projectors, one obtains, 

P^'^{r)P'''''{p)q''V^^,{q,r,p) = [m\r) - m\p)]P^\r)P'^-\p) . (3.15) 

On the other hand, contracting the full expansion of the vertex (13 ■2p by these transverse 
projectors and then contracting the result with the momentum of the background leg, we 
get 

q''P^'^{r)P^'''{p)V^,M. r,p) = iimSig,. + B^q^q,]P>^' ^{r)P^'^ [p) , (3.16) 
where the relation of Eq (13.111) has been used. Thus, equating both results, one arrives at 

iI{q)Bi{q,r,p) = m^{r) — m^{p) , B2{q,r,p) = 0. (3-17) 

The above relations, together with those of Eq. (13. lip , determine exactly the form factors 
Vi and V2 of the vertex Va^^, namely 

Vi(g,r,p) = '"'^''^~^'^^\ V2{q,r,p)=0. (3.18) 

We will now carry out the Taylor expansion of both sides of Eq (I3.17P in the limit g — )■ 0. 
To that end, let consider the Taylor expansion of a function f{q,r,p) around g = (and 
r = —p). In general we have 

f{q, -P -q,p) = f{-p,p) + [2(g ■ p) + qVi-p,p) + 2(g • P?n-P,p) + 0{q^) , (3.19) 

where the prime denotes differentiation with respect to {p + qY and subsequently taking the 
limit g — )■ 0, i.e. 



f(-p »)- limP^^^'~^~^'^U 
J[ p,p)- nm| g^p^gy j 



(3.20) 



Now, if the function is antisymmetric under p -H- r, as happens with the form factors -81,2, 
then f{—p,p) = 0; thus, for the case of the form factors in question, the Taylor expansion 
is {i = 1,2) 

Biiq, -p -q,p) = [2(g • p) + q^]B[{-p,p) + 2(g • pfB'l^-p^p) + ^(g^) . (3.21) 

Using Eq (I3.2ip . and the corresponding expansion for the rhs, 

m^{r)-m^{p) = m^{q+p)-m^{p) = 2{q-p)[m^{p)] + 0{q^) , (3.22) 

assuming that the /(O) is finite, and equating the coefficients in front of {q ■ p), we arrive at 
(Minkowski space) 

[m'{p)]' = tI{Q)B[{p). (3.23) 
13 



We emphasize that this is an exact relation, whose derivation rehes only on the WI and Bose- 
symmetry that Va^u{q,r,p) satisfies, as captured by Eq. fl2.14p and Eq. fl3.13p . respectively. 
The Euclidean version of Eq. (I3.23p is given in Eq. flS.ip . 

C. "One-loop dressed" approximation for the transition function 

We will next approximate the transition amplitude /^(g), connecting the gluon with the 
massless excitation, by considering only diagram (di) in Fig. IH corresponding to the gluonic 
"one- loop dressed" approximation; we will denote the resulting expression by 7a (g). 

In the Landau gauge, laiq) is given by 

Uq) = 1:Ca [ A(fc) A(fc + q)T^pxP^f'{k)P^''{k + q)B,,{-q, -k, k + q) , (3.24) 
^ Jk 

where the origin of the factor 1/2 is combinatoric, and Tai3x is the standard three-gluon 
vertex at tree-level, 

ra/.^(g, r,p) = g^^ij- - p)a + gau{p " q) ^, + QaM " r)^ . (3.25) 

To determine the corresponding transition function from Eq. f l3.10p . use that 

g"r,;3A(g, -k -q,k) = -{k + qy]g/sx + [{k + q)p{k + q)x - kpk^] , (3.26) 

to write 

I{q) = l^[k' -{k + q)']A{k)A{q + k)P^{k)P^%k + q) B,,{-q, k + q, -k) . (3.27) 

In the last step we have used the property of Eq. fl3.12p in order to interchange the arguments 
of B^^, so that the Taylor expansion of Eq. fl3.2ip may be applied directly; this accounts for 
the additional minus sign. Then, after the shift k + q ^ k, and further use of Eq. (13. 120 . 
/(g) becomes 

Hq) = / k'A{k)A{k + q)P^{k)P^%k + q) [B^g^, + B^q^^q,] . (3.28) 

1 Jk 

To obtain the limit of /(g) as g^ — )■ 0, we will employ Eq. f l3.2ip for Bi and /?2, as well as 
A{k + g) = A{k) + [2(g ■ k) + q^]A\k) + 2(g ■ kfA"{k) + ^(g^) . (3.29) 
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Observe that only the zeroth order term of P^u{k + q), namely Pfj_u{k), contributes in this 
expansion. Then, using spherical coordinates to write {q ■ k)^ = q^k"^ cos^ 9, and the integral 

/ f{k)cos^e = - [ f{k), (3.30) 

Jk " Jfc 

the /(g) in Eq. fl3.28p becomes in the limit — )■ (in d = 4) 

7(0) = -3CA^jk'A\k)B[{k) + \jk'^,[A\k)B[{k)]^ . (3.31) 

Then, partial integration yields 

jk^±-^[A\k)B[{k)] = -2 jk^A\k)B[{k), (3.32) 

and finally one arrives at (Minkowski space) 

m = \cAjk'A\k)B[{k). (3.33) 

The Euclidean version of this equation, Eq. (15. 3p . will be used in Section IVl 

We end this subsection with a comment on the dimensionality of the various form factors. 
The vertex Vaij.v has dimension [m], and so Vi, V2 and V3 are dimensionless, while the 
remaining form factors have dimension [m]~^. The integral /(g) has the same dimension as 
/?!, and as a result, in order to keep Vi dimensionless, Bi must have dimensions of [m]. 



D. Relating the gluon mass to the transition function 

In this subsection we show how the vertex V gives rise to a gluon mass when inserted 
into the corresponding SDE. We will restrict ourselves to the two diagrams shown in Fig. [U 
and will finally express m^(0) exclusively in terms of /(O), which, in turn, depends on the 
existence of Bfj_y through Eq. (13.331) . 

In the PT-BFM scheme, the SDE of the gluon propagator in the Landau gauge assumes 
the form 

A-'tf )/^.(.) ^ . (3.34) 

The most straightforward way to relate the gluon mass to the transition function / is to 
identify, on both sides of fl3.34p . the co-factors of the tensorial structure q^q^/q^ which survive 
the limit g^ — )■ 0, and then set them equal to each other. Making the usual identification (in 
Minkowski space) A~^(0) = — m^(0), it is clear that Ihs of f l3.34p furnishes simply 

\h.s\i^i. = m2(0) . (3.35) 
15 





FIG. 5: The "squared" diagram. 



It is relatively straightforward to recognize that the analogous contribution from the rhs 
comes from the standard "squared" diagram, shown in Fig. [51 Specifically, the starting 
expression is 

n'^'^(g) = Ig'CA I T^^^P'''{k)P^^{k + q)[T + V];^A{k)A{k + g) + ■ ■ ■ , (3.36) 

where, as explained earlier, the (all order) vertex JT has been replaced by F + and the 
ellipses denote terms that, in the kinematic limit considered, do not contribute to the specific 
structure of interest. 

The relevant contribution originates from the part containing the vertex V, to be denoted 
by Il^'^{q)\y; it is represented by the diagram in Fig. [5l In particular, by virtue of Eq. (13. 6p . 
we have 

n^'(g)|v = ^g'CAl^A{q + k)A{k)T^^^P"''{k)P^^{k + q)U;^ 

= 9'{^lcAj^A{q + k)A{k)T^^^P-^{k)P^'^{k + q)B,^^ (^) ^^^^^ 

= ^^-^9'I^g), (3.37) 
q^ 

where in the second line we have used Eq. (13.71) [ with Iu{q) — ?• Iu{(l)], and Eq. (13.241) in the 
third line. 

n 

Thus, using the fact that, since L(0) = [42], from the identity of Eq. (I2.13p we have 
that 1 + G(0) = F~\0), then the rhs of iKM} becomes 

rhs|.M9. = -/F2(0)P(0). (3.38) 
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We next go to Euclidean space, following the usual rules, and noticing that, due to the change 
= i J^^ we have P{0) — t- — /|(0); so, equating (I3.35P and (I3.38P we obtain (suppressing 
the index "E") 

m\0) = g^F\0)P{0) . (3.39) 
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Note that the m^(0) so obtained is positive-definite. We emphasize that the relation of 
Eq. f l3.39p constitutes the (gluonic) "one-loop dressed" approximation of the complete rela- 
tion; indeed, both the SDE used as starting point as well as the expression for / are precisely 
the corresponding "one-loop dressed" contributions, containing gluons (but not ghosts). 
Finally, let us consider the exact relation 15 1 

m\q^) = \l + G{q^)fm\q^), (3.40) 

expressing the dynamical mass m{q^) of the standard gluon propagator A(g^) in terms of the 
corresponding mass, m{q^), of the PT-BFM gluon propagator [usually denoted by A(g^)] 
in the same gauge [in this case, in the Landau gauge]. At = this relation reduces to 
m^(0) = m^(0)F^(0), so that Eq. (I3.39P may be alternatively written as 

m\<S) = g^P{0) . (3.41) 

Interestingly enough, when written in this form, the mass formula derived from our SDE 
analysis coincides with the one obtained for the photon mass in the Abelian model of Jackiw 
and Johnson (Eq. (2.12) in [isl). In addition, this last form facilitates the demonstration of 
the decoupling of the massless excitation from the on-shell four-gluon amplitude (see Section 

lYiD. 

In principle, the analysis presented above may be extended to include the rest of the 
graphs contributing to the gluon SDE, invoking the corresponding pole parts of the remaining 
vertices; however, this lies beyond the scope of the present work. 



IV. BS EQUATION FOR THE BOUND-STATE WAVE FUNCTION 

As has become clear in the previous section, the gauge boson (gluon) mass is inextricably 
connected to the existence of the quantity B[. Indeed, if B[ were to vanish, then, by virtue of 
f l3.33p so would /(O), and therefore, through fl3.39p we would obtain a vanishing m^(0). Thus, 
the existence of B[ is of paramount importance for the mass generation mechanism envisaged 
here; essentially, the question boils down to whether or not the dynamical formation of a 
massless bound-state excitation of the type postulated above is possible. As is well-known, 
in order to establish the existence of such a bound state one must (i) derive the appropriate 
BSE for the corresponding bound-state wave function, S^,^, (or, in this case, its derivative), 
and (ii) find non-trivial solutions for this integral equation. 

17 



{«) (6) 




FIG. 6: The complete BSE for the fuh three gluon vertex Jr'^^j^{q,r,p). 

To be sure, this dynamical equation will be derived under certain simplifying assumptions, 
which will be further refined in order to obtain numerical solutions. We emphasize, therefore, 
that the analysis presented here is meant to provide preliminary quantitative evidence for the 
realization of the dynamical scenario considered, but cannot be considered as a conclusive 
demonstration. 

The starting point is the BSE for the vertex T'^^j^{q,r,p), shown in Fig. [61 Note that, 
unlike the corresponding SDK of Fig. [21 the vertices where the background gluon is enter- 
ing (carrying momentum q) are now fully dressed. As a consequence, the corresponding 
multiparticle kernels appearing in Fig. [6l are different from those of the SDE, as shown in 
Fig. [3 




FIG. 7: (A) Schematic relation between the SDE and BSE kernels. (B) Example of a diagram 
not contamed in the corresponding BSE kernel, in order to avoid over counting. 
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(A) cWW^• = cWWX^ + cWWN^ 






m, n ► m , 

k + a 



FIG. 8: (A) The separation of the vertex in regular and pole parts. (B) The BSE for the bound- 



state wave function B^y. 



The general methodology of how to isolate from the BSE shown in Fig, pthe correspond 
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22|. Specifically, 



ing dynamical equation for the quantity B^^ has been explained in 
one separates on both sides of the BSE equation each vertex (black circle) into two parts, a 
"regular" part and another containing a pole 1/g^; this separation is shown schematically in 
Fig. [HI Then, the BSE for B^,^{q, r,p) is obtained simply by equating the pole parts on both 
sides. Of course, for the case we consider the full implementation of this general procedure 
would lead to a very complicated structure, because, in principle, all fully dressed vertices 
appearing on the rhs of Fig. [6] may contain pole parts [i.e., not just the three-gluon vertex 
of (a) but also those in (b), (c), and (d)]. Thus, one would be led to an equation, whose Ihs 
would consist of B^,^, but whose rhs would contain the B^^^, together with all other similar 
vertices, denoted by in Eq. (13. 14^ . Therefore, this equation must be supplemented by 
a set of analogous equations, obtained from the BSEs of all other vertices appearing on the 
rhs of Fig. |6] [i.e., those in (b), (c), (d) ]. So, if all vertices involved contain a pole part, one 
would arrive at a system of several coupled integral equations, containing complicated com- 
binations of the numerous form factors composing these vertices (see, for example. Fig. 11 



m 



221). 

It is clear that for practical purposes the above procedure must be simplified to something 
more manageable. To that end, we will only consider graph (a) on the rhs of Fig. |6l thus 
reducing the problem to the treatment of a single integral equation. 
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^'P '".M (qi) (02) (as) 



FIG. 9: The Feynman diagrams considered for the BS kernel. The interaction vertices are approx- 
imated by their tree level values, while the internal gluon propagators are fully dressed. 

Specifically, the BSE for B^^, is given by [see Fig. E] 

Jk 

In addition, we will approximate the four-gluon BS kernel K, by the lowest-order set of 
diagrams shown in Fig. [9], where the vertices are bare, while the internal gluon propagators 
are fully dressed. 

To proceed further, observe that the diagram (ai) does not contribute to the BSE, because 
the color structure of the tree-level four-gluon vertex vanishes when contracted with the 
color factor /"'"^ of the Ba/3. Diagrams (02) and (03) are equal, and are multiplied by a Bose 
symmetry factor of 1/2. So, in this approximation, the BS kernel is given by 

/C~(-fc,p, -p-q,k + q) = -ig'r^n^T^^lA^\k - p)T'^X ' (4-2) 

where F*^*^^ is the tree-level value of the three gluon vertex. So, using this kernel and setting 
the gluon propagators in the Landau gauge, the BSE becomes 

B^, = -2TiiasCA [ B^pA{k + q)A{k)A{k-p)P''''{k + q)P^''{k)P'^\k~p)T^^lT^';l (4.3) 
Jk 

where we have cancelled out a color factor /"'"^ from both sides. 

Let us focus on the Ihs of Eq. fl4.3p . Using the Taylor expansion in Eq. (13.211) . the fact 
that B2 = [see Eq. (I3.17p ]. and multiplying by a transverse projector we obtain, 

P^''{p)B^, = 6{q-p)B[{p) + 0{q') . (4.4) 

Next, let us denote by [rhs]^,y the rhs of Eq. (14.31) . Inserting the bare value for the three 
gluon vertices, multiplying by the transverse projector, and using the Taylor expansions in 
Eq. (I3.2ip and (13.290 . after standard manipulations one obtains the result 

P>'''{p)[ThsU = -4ma,CA{q-p) [ B[{k)A\k)A{k - p)Af{p, k) + 0{q^) , (4.5) 

Jk 
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where we have defined the kernel 



[8p'e - 6{pk) ip' + A;') + 3(/ + r ) + (pkY] . (4.6) 



p4fc2(A;-p)2 

Thus, equating the Ihs with the rhs, we derive the following BSE for the derivative of the 
form factor that appears in the mass relation Eq. ( 13.23p . 



B[{p) 



2Tci 



^ asCA / B[{k)/\\k)/\{k-p)M{p,k) 
-J Jk 



Going to Euclidean space, we define 

X = p^ ] y = k'^ ; z = {p + k^ , 
and write the Euclidean integration measure in spherical coordinates. 



d'^kE 1 



(27r)4 {2ny 



dyy / d9 sin^ 9 . 



(4.7) 



(4.8) 



(4.9) 



so that the BSE becomes 

B[{x) = / dyyB[{y)A\y) 



d6 sin 6 cos 6 



z + 10(a; + y) + -(a;^ + + lOxy) 



In spherical coordinates we have that z = x + ?/ + 2^fxycosQ. So, around x = 0, 

1 1 



z X + y 



jxy 
l-2^^cos^ 

X + y 



(4.10) 



(4.11) 



and using the Taylor expansion for the gluon propagator A{z), the limit x — > can be taken 
in the BSE, giving the value 



B[{0) = lim5;(x) 



asCA 



dyy'B[{y)A\y)A'{y). 



(4.12) 



x->0 ' StT Jo 

Let us finally implement an additional simplification to Eq. fl4.10p . which will allow us 
to carry out the angular integration exactly, thus reducing the problem to the solution of a 
one- dimensional integral equation. Specifically, the simplification consists in approximating 
the gluon propagator A{z) appearing in the BSE of fl4.10p [but not the A^(?/)] by its tree 
level value, that is, A{z) = 1/z. Then, with the aid of the angular integrals, 
, sin^ 9 cos 9 



d9- 



ly I " ^sm'^9cos9 
- / d9 



TT 

16x 

TT 

Ax 



y 

— {y - 2x)Q{x -y) + -^{x - 2y)Q{y - x) 
x^ y^ 



y *^ 

— e(x -y) + -^<S>{y - x) 

x^ 2/ 



(4.13) 
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one brings Eq. fl4.10p into the form 

2hy 3y2 



Ax Ax 



+ 



+ jyyB[(y)A\y)v{3 + ^''--\fjy (4.14) 

The hmit a; ^ of this equation is given by (the change of variable y = tx may be found 
useful) , 

POO 

B[{0) = ^ / dyyB[{y)A\y) . (4.15) 

OTT Jo 

Note that this result coincides, as it should, with that obtained from Eq. (I4.12p after setting 
^'{y) = ^V?/^5 namely the derivative of the tree-level propagator. 



V. NUMERICAL SOLUTIONS AND EXISTENCE OF A BOUND-STATE 

In this section we will carry out a detailed numerical analysis of the integral equation 
obtained in the previous section, namely Eq. fl4.14p . 

First of all, let us point out that, despite appearances, the integral equation fl4.14p is 
not linear in the unknown function B[{x). The non-linearity enters through the propagator 
A{y), which depends on the dynamical mass m^{y) through Eq. (12. 6p : as a result, and by 
virtue of Eq. (13.231) . which, in Euclidean space reads 

[m'iy)r = -IiO)B[iy), (5.1) 

it is clear that A{y) depends on B[{x) in a complicated way. Specifically, from the two 
aforementioned equations we have 



A-\y) = yJ{y) + m^{y), 



y 



m\y) = m\0)-I{0) / dzB[{z) . (5.2) 



where J(0) may be approximated by its "one-loop dressed" version /(O) given in (13.331) . 
which in Euclidean space becomes 

o/^ roc 

/(0) = |^ j^dyy'A\y)B[{y). (5.3) 

Then, Eq. (I4.14p must be solved together with the two additional relations given in 
Eq. (15. 2p . as a non-linear system. In addition, one may use Eq. (13.390 . in order to obtain an 
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(approximate) constraint for /(O). Note also that Eq. fl4.14p . again due to Eq. (15. ip . may be 
recast entirely in terms of m'^{y) and its derivative. 

For the purposes of the present work we will simplify somewhat the procedure described 
above. Specifically, we will present two different approaches, each one particularly suited for 
probing distinct features of Eq. (I4.14p and the accompanying Eqs. (15. 2p . In particular, we will 
first study Eq. (I4.14p in isolation, using simple Ansatze for A{y). The purpose of this study is 
to establish the existence of non-trivial solutions for B[, study their dependence on the value 
of the strong coupling as, and verify the asymptotic behavior predicted by Eq. (15. 4p . Of 
course, since A(?/) at this level is treated as an "external" object, the homogeneous Eq. (I4.14p 
becomes linear in B[; as a result, given one solution we obtain a family of such solutions, 
through multiplication by any real constant. Then, as a second step, we will use the available 
lattice data for the gluon propagator A(?/), in order to obtain the corresponding solution for 
B[. Now, the linearity induced by treating A{y) as an external input will be resolved by 
resorting to Eq. (15. 2 p and Eq. ( ]3.39p : thus, out of the infinite family of solutions only one 
will be dynamically selected. These two approaches will be presented in subsections IV Bl 
and IV CI while subsection IV Al deals with the asymptotic behavior of the solutions. 



A. Asymptotic behavior 

Before turning to the numerical treatment of Eq. (I4.14p . it is useful to study the behavior 

of the solutions for asymptotically large values of x. Setting A{y) = 1/y, it is relatively 

straightforward to establish that the equation admits a power-law solution of the form 

B[{x) = Ax''. Specifically, substituting this Ansatz into the first integral of Eq. (I4.14p . 

which is the dominant part for large x, and carrying out the integrations, one arrives at the 

following algebraic equation for b, 

2Att _ 3 25 3 

~ b + i ^ 4(6 + 2) " 4(6 + 3) ' 

together with the restriction 6 > — 1, imposed in order to assure convergence in the lower 
[y = 0) limit of integration. Setting A = 2A7T/asCA, one arrives at the third-order equation 

4A6^ - (34 - 24A)6^ - (151 - 44A)6 - 141 + 24A = , (5.5) 

that may be easily solved; the solution that satisfies 6 > — 1 is shown in Fig. [10] as a function 

of ttg. 
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FIG. 10: The physically relevant solution of Eq. (15. Sp . 
B. The linearized case: solutions for various gluon propagators 



Next we discuss the numerical solutions for Eq. fl4.14p for arbitrary values of x. Evidently, 
the main ingredient entering into its kernel is the nonperturbative gluon propagator, A(g). 
In order to explore the sensitivity of the solutions on the details of A(g), we will employ 
three infrared- finite forms, to be denoted by Ai(g), A2(g), and A3(g), focusing on their 
differences in the intermediate and asymptotic regions of momenta. 

Let us start with the simplest such propagator, namely a tree-level massive propagator 
of the form 

A^\q^) = q^ + ml, (5.6) 

where mg is a hard mass, that will be treated as a free parameter. On the left panel of 
Fig. [m the (blue) dotted curve represents Ai(g^) for mo = 376 MeV. 

The second model is an improved version of the first, where we introduce the 
renormalization-group logarithm next to the momentum g^, more specifically 



1 + 



In 



+ pml 



+ ml. 



(5.7) 



967r^ V 

where p is an adjustable parameter varying in the range of p G [2, 10]. Notice that the hard 
mass appearing in the argument of the perturbative logarithm acts as an infrared cutoff; 
so, instead of the logarithm diverging at the Landau pole, it saturates at a finite value. The 
(black) dashed line represents the Eq. (15. 7p when p = 16, ttiq = 376 MeV, and p = 4.3 GeV. 
The third model is simply a physically motivated fit for the gluon propagator determined 
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FIG. 11: Left panel: The three models for the gluon propagator as function of the momentum 
g^. The (red) continuous hne is the fit for the lattice gluon propagator given by Eq. ()5.8p when 
m = 520 MeV, = 5.68, pi = 8.55, p2 = 1.91, and p = 4.3 GeV; the (black) dashed line is the 
model of Eq. (j5.7p with p = 16, Os = 0.667 and tjiq = 376 MeV, while the (blue) dotted line rep- 
resents the massive propagator of Eq. (15. 6p when mo = 376 MeV. Right panel: The corresponding 
solutions of Eq. ()4.14p obtained with the gluon propagators shown on the left panel. The solutions 
for B[{q) are obtained when we fix the value of ct^ = 1.48, = 0.667, and as = 0.492 for Ai(g), 
A2{q), and A3 ((7), respectively. 

by the large- volume lattice simulations of Ref . \q] , and shown on the left panel of Fig. [HI 
The lattice data presented there correspond to a SU{3) quenched lattice simulation, where 
A(g) is renormalized at /i = 4.3 GeV. This gluon propagator can be accurately fitted by the 
expression (e.g.,{l^) 



where m^iq^) is a running mass given by 



(5.8) 
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^M) = ^ ^, (5.9) 

and the values of the fitting parameters are m = 520 MeV, (7^ = 5.68, pi = 8.55 and, 
p2 = 1.91. On the left panel of Fig. [TTl the (red) continuous line represents the fit for 
the lattice gluon propagator given by Eq. (15. 8p . Notice that, in all three cases, we have fixed 
the value of A^^(O) = ml ^ 0.14. 
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FIG. 12: The (black) continuous curve represents the solution obtained from Eq. ()4.14p using the 
propagator Ai{q) of Eq. (|5.6p . with a<j = 1.48 and mo = 376 MeV. The (red) dotted line is the 
best fit obtained for the asymptotic behavior of B[{x) given by B[{q) = Aq^^ with Ai = 14.80 and 
b = —0.756. Notice that this value is in excellent agreement with the power found by the analytical 
determination, shown in Fig. [lOl 

Our main findings may be summarized as follows. 

Fig. [m right panel, we show the solutions of Eq. (I4.14p obtained using as input the 
three propagators shown on the left panel. For the simple massive propagator of Eq. (15. 6p . a 
solution for B[{q) is found for = 1.48; in the case of A2(g) given by Eq. (15. 7p . a solution 
is obtained when = 0.667, while for the lattice propagator A3(g) of Eq. (15. 8p a non-trivial 
solution is found when as = 0.492. 

(ii) Note that, due to the fact that Eq. (I4.14p is homogeneous and (effectively) linear, 
if B[{q) is a solution then the function cB[{q) is also a solution, for any real constant c. 
Therefore, the solutions shown on the right panel of Fig. [11] corresponds to a representative 
case of a family of possible solutions, where the constant c was chosen such that B[{0) = 1. 

(Hi) Another interesting feature of the solutions of Eq. (I4.14p is the dependence of the 
observed peak on the support of the gluon propagator in the intermediate region of momenta. 
Specifically, an increase of the support of the gluon propagator in the approximate range 
(0.3-1) GeV results in a more pronounced peak in B[{q). 

(iv) In addition, observe that due to the presence of the perturbative logarithm in the 
expression for A2(q') and A3(g), the corresponding solutions B[{q) fall off in the ultraviolet 
region much faster than those obtained using the simple Ai(g) of Eq. (15.61) . In order to 
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FIG. 13: Left panel: The behavior of the gluon propagator Ai{q), given by Eq. (I5.6p . for various 
values of rriQ in the range of 300 — 800 MeV. Right panel: The corresponding solutions for B[{q), 
obtained using the gluon propagators shown on the left panel. For each value of niQ, we found that 
the solution for B[ (g) is obtained for a particular value of ag ■ 

check whether the power-law asymptotic behavior, B[{q) = Aq"^^, determined in our previous 
analysis, is in agreement with our numerical solution, we isolate in Fig. [12] the solution of 
B[{q) obtained with Ai(g) and = 1.48 (black continuous curve) and compare it with the 
best fit obtained for large values of q^ (red dotted curve). Indeed, the asymptotic tail of 
B[{q) falls off as power law of the type B[{q) = Aq^^ with A = 14.80 and 6 = -0.756. Notice 
that the value of b obtained from the fit is in perfect agreement with values obtained from 
Eq. (15. 5p , shown in Fig. [101 

(v) On the left panel of Fig. [13] we plot Ai(g), given by Eq. (15. 6p . for different values of 
mo in the range of 300 — 800 MeV. In order to determine how the solutions are modified 
when one varies the value of mo, we show on the right panel of Fig. [13] the various B[{q), all 
of them normalized at -8^(0) = 1. As we can see, the solutions display the same qualitative 
behavior; however, for each mo, the non-trivial solution is obtained for a different value of 
as- In fact, as the values of mo increase, so do the values of as needed for obtaining a 
solution; the exact dependence of as on ml is shown in Fig. [TU 

(vi) Next, we study how size variations in the intermediate region of the gluon propagator 
change the values of ag needed in order to obtain non-trivial solutions from Eq. (I4.14p . To 
address this point systematically, we employ the gluon propagator A2(g) of Eq. (15. 7p . varying 
the parameter p in the range of p G [2, 16], keeping fixed mo = 367 MeV, as shown on the left 
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FIG. 14: The values of Og furnishing non-trivial solutions to Eq. (|4.14p as we vary tjiq in the Ai((7) 

of Eq. (fsrei) . 

panel of Fig. [151 the corresponding B[{q) for each value of p are plotted on the right panel. 
Evidently, decreasing p increases the support of the gluon propagator in the intermediate 
region, and, as a result, one needs smaller value of in order to obtain solutions for B[{q). 
This last property is better seen in Fig. [161 where we present the values of needed to 
solve Eq. f l4.14p as one varies p in A2(g). 



2.0 




I 1 — I 1 — I 1 — I 1 — I I ' — I 1 — I ' — I ' — I 1 — I 

1E-3 0.01 0.1 1 10 1E-3 0.01 0.1 1 10 100 

q'[GeV'] q'[GeV'] 



FIG. 15: Left panel: The behavior of the gluon propagator, A2{q), given by Eq. (|5.7p . when the 
value of mo = 367 MeV is fixed, and p varies in the range 2-16. Right panel: The corresponding 
solutions for B[{q) obtained with the gluon propagators shown on the left panel. To each value of 
p corresponds a specific value of that yields a solution B[{q). 



28 



0.65- 
0.60- 




































■ 


0.55- 
0.50- 





































































2 4 6 



10 12 14 16 18 



FIG. 16: The values of Og for which we obtain non-trivial solutions to Eq. ()4.14p as we vary p in 
A2{q) of Eq. 

C. Non-linear treatment: uniqueness of B[{x) and 



In the previous subsection we have practically solved Eq. f l4.14p in isolation, in the sense 
that we have not used the supplementary conditions of Eq. (15. 2p . and have treated A(g) as 
an external independent quantity. As a result, the homogeneous Eq. (14. 14p was effectively 
linearized, giving rise to families of solutions cB[{x), parametrized by the value of c. In this 
subsection we will restore the non-linearity of Eq. (I4.14p : as a result, the arbitrariness in 
the value of c is completely eliminated, and one obtains a single expression for B[{x) and 
m^(x), for a unique value of 

The way a unique solution for B[{x) is singled out [i.e., a value for c is dynamically 
chosen] is by combining Eq. (I3.39P and Eq. (15. 3p : specifically, we will require that the value 
for 7(0) obtained from the former equation coincides with that obtained from the latter, 
namely that 

^r. r°° _ „ 

(5.10) 



v/A-i(0)/4vra,F2(0) = |^ / dy A' (y) B[{y) 



Now, the Ihs of (I5.10p is fixed, because, as mentioned in the previous subsection, we must 
have as = 0.492 in order for Eq. (I4.14p to have solutions for this particular (lattice) prop- 
agator as input, while A~^(0) and F{0) are fixed from the lattice. Specifically, the SU{3) 
large-volume lattice simulations of Ref. ^| yield ^ 0.141 (see Fig. [TT]) and F{0) ~ 2.76 
(see Fig. [T71) . Then, the integral on the rhs (I5.10p must match the value of the Ihs, and this 
can only happen for one particular member of the family cB[{x). 

In Fig. [181 we show the solution for B[[x), which satisfies the constraints imposed on the 
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Ghost dressing function F(q} 
n = 4.3 GeV 
Lattice p = 5.7 

is L=64 
+ L=SO 
Fit 
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n 

FIG. 17: Lattice results |6i] for the ghost dressing function, F(q), renormalized at /x = 4.3 GeV. 
Notice that F(0) w 2.76. 



value of /(O), obtained when a, = 0.492 and B[{0) = 0.086. 
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FIG. 18: The numerical solution B[{q) obtained from Eq. (j4.14p . under the constraints imposed 
by Eqs. (ICTl) and ([3331), and with = 0.492. 



Once the unique solution B[{q) has been determined, one may use Eq. (15.20 to determine 
the behavior of the squared gluon mass m^(x). Integrating numerically B[[q) and fixing 
m(0) = 0.14, we obtain the result shown in Fig. [191 

Evidently, the function rn?{q^) displays a plateau in the deep infrared, and then decreases 
sufficiently fast in the ultraviolet region, as expected on general grounds 



30 



> 

C3 



1 1 1 1 1 1 1 

Dynamical gluon mass 



























































FIG. 19: The square of the dynamical gluon mass, obtained from Eq. (j3.23p . after plugging into it 
the B[ (x) shown in Fig. [TSl 

VI. DECOUPLING OF THE MASSLESS EXCITATION: AN EXAMPLE 

In this section we give an explicit example of how the massless excitation decouple from 
an on-shell amplitude. Specifically, we will show how this is indeed what happens in the 
case of the four-gluon amplitude. To be sure, a complete proof of the decoupling of the 
massless excitation from all Yang-Mills amplitudes requires the treatment of kernels with 
an arbitrary number of incoming gluons, which is beyond our powers at present. However, 
the example considered here captures the essence of the underlying decoupling mechanism. 



The demonstration followed here is similar to that given in [19] for the case of an Abelian 
model. One starts by considering the complete four-gluon amplitude, [graph (a) in Fig. 120] . 
which consists of three distinct pieces: (i) the amplitude represented by the diagram (&), 






ib) (c) id) 

FIG. 20: The complete four-gluon amplitude and the various terms composing it. 
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which is regular as — )■ 0, (ii) the graph (c), which contains the massless excitation, 
coupled to the external gluons through the proper vertex function 5, and (iii) the one- 
particle reducible term, denoted by (rf), which is excluded from the SDE kernel in the usual 
skeleton expansion. Of course, the above amplitudes are none other than (62), (^3), and 
(a) in Fig. [31 respectively. Since the amplitude (6) is regular by construction, one must 
only demonstrate that, as — )■ 0, the divergent part of (c), whose origin is the massless 
excitation, cancels exactly against an analogous contribution contained in (ci), leaving finally 
a regular result. 

We start by considering the term {d). Within the PT-BFM framework that we use, 
the off-shell gluon (carrying momentum q) is effectively converted into a background gluon; 
thus, the gluon propagator appearing inside {d) is given by A(g), while the two three-gluon 
vertices are the F' defined in Eq. (12. 9p . So, 

{d) = -ig^r^^M.pi.P2)P''^{q)Hq)Kp^{q.p^.PA) 

= -^/r^/.z.(?,Pi,P2)A(g)r'^^(g,P3,P4) , (6.1) 

where the factor {—i) comes from the definition of the gluon propagator, Eq. (12. ip . In 
the second line we have eliminated the longitudinal term inside P'^^{q) using the 

"on-shellness" condition 

= 0, (6.2) 

valid for both three-gluon vertices. We emphasize that the full F' is needed (with the V 
included) for the on-shellness condition of Eq. (16. 2 p to be fulfilled. Note also that, if one 
had chosen a non vanishing gauge-fixing parameter ^ for the gluon propagator (instead of 
the ^ = of the Landau gauge), then the condition of Eq. (16. 2p is instrumental for the 
cancellation of the unphysical parameter ^ from the physical amplitude. 

Next, it is clear that from the vertex V contained in F' only the part U survives, [see 
Eq. (13. 3p ]. because all longitudinal momenta contained in R are annihilated on shell, i.e. 
when contracted with the appropriate polarization vectors e^{p), due to the validity of the 
relation p^e^(p) = 0. Then, we have that (suppressing indices) 

r'Ar' = (r + f/)A(r + f/) 

= TAT + T'AU + UAT' -UKu . (6.3) 
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Given that the first term in (16.31) is regular, while the second and third term vanish on shell 
by virtue of (16. 2 p [which is triggered because Uafiu is proportional to qa, see Eq. (13. 4p ]. we 
are led to the following expression for the pole part of (d) 

{d)poie = ig'Uo.f.uA{q)U^^. (6.4) 

Then, using Eqs. (13. 7p and (13. 9p . we obtain 

{dUo = - {i? i?} [//^(g)A(g)] . (6.5) 

Now, in the limit — )■ 0, the quantity in square brackets goes to 1, precisely by virtue of 
Eq. (I3.4ip [remember, A~^(0) = m^(0)]. Therefore, 

lim(dW = -lim{B(^)B}. (6.6) 



which is precisely the contribution of the term (c) in the same kinematic limit, but with the 
opposite sign. Therefore, the on-shell four-gluon amplitude is free from poles at = 0, as 
announced. 

Finally, note that, as an alternative, one might opt for eliminating completely any refer- 
ence to V in the amplitude {d) from the very beginning, namely the first step in Eq. (16.11) : 
this is of course possible, given that some parts of the fully longitudinal vertex V vanish on 
shell, while the rest vanishes when contracted with the transverse projector Pap{q). In such 
a case, however, one may not dispose of the longitudinal part of Papiq) any longer, because 
now the on-shellness condition of Eq. (16.21) is distorted, precisely due to the absence of V . It 
is a straightforward exercise to demonstrate that if one were to take the produced mismatch 
into account, one would recover exactly the same result found above. 



VII. DISCUSSION AND CONCLUSIONS 



The gauge-invariant generation of a gluon mass relies on the existence of massless bound- 
state excitations, which trigger the Schwinger mechanism. The presence of these excitations 
in the skeleton expansion of the full three-gluon vertex F^^,^ induces longitudinally coupled 
pole structures, giving rise to a purely nonperturbative component, the pole vertex Vaiiv 

In this article we have studied in detail the dynamical ingredients associated with the 
vertex V; in particular, the poles in V are identified with the propagator of the massless 
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scalar excitation, while the tensorial structure is determined by two basic purely nonpertur- 
bative quantities: the transition amplitude, denoted by which at the diagrammatic level 
connects the gluon propagator with the massless scalar propagator, and the effective vertex, 
denoted by B^^, connecting the massless excitation to two outgoing gluons. 

The powerful requirement of maintaining the gauge invariance of the theory intact re- 
stricts the form of the various form factors composing B/j^^, and establishes a non-trivial con- 
nection between the transition function and the first derivative of the momentum-dependent 
gluon mass. The insertion of the vertex Va^u (or, effectively, its surviving component Uafi^) 
into the SDE of the gluon propagator (in the Landau gauge) allows one to express, at zero 
momentum transfer, the gluon mass in terms of the transition function, by means of a rather 
simple formula. In fact, it turns out that, the relevant dynamical quantity is the derivative 
of the form factor Bi, denoted by B[. 

As we have demonstrated, in the aforementioned kinematic limit, the homogeneous BSE 
obeyed by the reduces in a natural way to an analogous integral equation for B[. The 
detailed numerical study of an approximate version of this latter equation reveals the exis- 
tence of non-trivial solutions for B[, which, when inserted into the corresponding formulas, 
furnish the momentum dependence of the gluon mass. The existence of these solutions adds 
weight to the hypothesis that the nonperturbative Yang-Mills dynamics lead indeed to the 
formation of the required massless bound-states. 

It is clear that some of the dynamical aspects of this problem merit a further detailed 
study, due to their relevance in the ongoing scrutiny of the infrared properties of the Yang- 
Mills Green's functions. Particularly important is to consider the effects of bound-state poles 
in the SD kernels of not only the three-gluon vertex, as we did in this article, but of all other 
fundamental vertices of the theory. Such an investigation would involve some or all of the 
vertices appearing in Eq. fl3.14p . which would form a coupled system of homogeneous integral 
equations. Given the recent lattice results on the ghost propagator, especially interesting in 
this context is the dynamical information that one might be able to obtain about the quantity 
B, corresponding to the wave-function of the ghost-ghost channel [(vertex Oi in Fig. (jl])]. 
Specifically, while the ghost dressing function F is found to be finite in the infrared, a fact 
that can be explained by the presence of the gluon mass that saturates the corresponding 
perturbative logarithm, there is no dynamical mass associated with the ghost field. One 
would expect, therefore, that the solution of the corresponding system should give rise to a 
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non- vanishing B[, as before, but to a vanishing ghost-ghost wave function B. 

As has been explained in detail in Section |Tll the incorporation the massless excitations 
modify the three-gluon vertex F, giving rise to the new vertex F', defined in Eq. fl2.9p . It 
would certainly be particularly interesting to compare the characteristic features of F' with 
results obtained on the lattice for the three-gluon vertex 43(]. In particular, one might in 
principle be able to relate the presence of the massless poles to possible divergences of some 
of the form factors, in the appropriate kinematic limit. To that end, one must first determine 
the closed expression for the entire vertex V from Eq. (12. 8p and the WI and STIs it satisfies. 



Then the answer should be written in a standard basis, such that of Ball and Chiu |44, |45| . 
and the final result projected on the particular kinematic configurations usually employed 
in the lattice calculations. We hope to be able to carry out this project in the near future. 
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